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We derive the gluon initiated contribution to the initial condition for the perturbative fragmen- 
tation function of a heavy quark through order 0{a 2 s ) in the MS scheme. This result is needed 
for the resummation with next-to-next-to-leading logarithmic accuracy of quasi-collinear logarithms 
ln fc (m 2 ) in heavy quark differential distributions by solving the complete DGLAP equation. Together 
with the previously evaluated fermion initiated components, this result completes the derivation of 
the initial condition for the perturbative fragmentation function at next-to-next-to-leading order. 

I. INTRODUCTION 

Production of heavy flavors (charm and bottom) in high energy processes has become an increasingly important 
subject in the last decade. In this type of processes, not-completcly inclusive observables like energy or transverse 
momentum distributions of heavy flavored hadrons are of special interest. The reason for that lies in their dependence 
on both the short-distance perturbative physics and the long-distance physics governing the formation of the heavy 
flavored hadrons. A peculiar feature of these observables is that they contain logarithmic sensitivity to the mass of 
the heavy flavor to all orders in perturbation theory. For that reason perturbative calculations of such cross-sections 
are not reliable when the typical hard scale in the process is much larger than the mass of the produced heavy quark. 

A formalism for the study of differential distributions for heavy quark production was proposed long ago by Mele 
and Nason . It is known as the Perturbative Fragmentation Function (PFF) formalism and its idea is that in order 
to restore the applicability of perturbative QCD to the type of observables discussed above, one has to resum classes 
of large logarithms to all orders in perturbation theory The method is process independent and is based on the 
DGLAP evolution equation 

The main goal of this paper is to derive the gluon initiated component of the initial condition for the perturbative 
fragmentation function at order 0{a 2 s ). Combined with the fermion initiated contributions Q at order 0{a 2 s ) 1 and 
with the three- loop time-like splitting functions (when the latter become available), this result can be used to perform 
resummations of the large collinear logs with next-to-next-to-lcading logarithmic (NNLL) accuracy. 

In general, the importance of the gluon initiated contribution increases in processes with larger separation between 
the hard scale and the quark mass or in processes where gluon production is not suppressed compared to the production 
of heavy quarks. This expectation is supported from the previous applications of the PFF formalism at the NLL 
accuracy level. The effect of the inclusion of the gluon component is certainly important for hadron colliders 0, || . 
The effect of the gluon component in e + e~ has been discussed, for example in [(|. Therefore, at the NNLL precision 
level, both the effect of the gluon initiated component and the quark-gluon mixing will be important. In fact, as we 
demonstrate in this paper, the property of the initial condition for the PFF - that its gluon initiated component at 
order 0(a s ) vanishes for fio = m, does not persist at order 0(a 2 ). 

This paper is organized as follows: in the next section we review some general results and then describe the 
derivation of the gluon initiated component. In section III we present our result and discuss its properties. At the 
end we present our conclusions. 

II. PROCESS-INDEPENDENT DERIVATION OF D'f 

Let us consider heavy flavor production in a hard scattering process. We are interested in the processes for which the 
characteristic hard scale Q is much larger than the mass m of the heavy quark. It follows from the QCD factorization 
theorems 0, H, 0, E3 that the non-perturbative contributions to the energy distribution of the heavy-flavored hadron 
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produced in such reaction are contained in the so-called non-perturbative fragmentation function. That function 
describes the formation of the observed heavy flavored hadron from the heavy quark Q at a scale set by the mass to, 
and is typically extracted from e + e~ data. In this paper we will not be concerned with the non-perturbative component 
of the fragmentation function; details about its implementation can be found in hlkl 5, 6, 11, 12, 13, 14, 15lhtf,hJ. 

We will be only interested in the perturbative production of a heavy quark with mass to << Q. In such kinematics 
all power corrections ~ (m/Q) p are neglected. This means that at large scales the heavy quark behaves as massless 
and the role of its mass is, essentially, to regulate the appearing collinear divergences. 

For the energy distribution of the heavy quark one writes: 

d<jQ . sr-^ f 1 dx da a 



dz 



Here to is the mass and z = Eq/ EQ^ max the energy fraction of the heavy quark produced in the reaction. The sum 
in Eq.JQ) runs over all partons - quarks, antiquarks (including the flavor Q) and gluons - that can be produced in 
the hard process and fj, is the factorization scale. The function da a /dx, also known as coefficient function, is the 
collinearly renormalized differential cross-section for producing a massless parton a. The coefficient function is related 
to the (bare) cross-section a a (z,Q,e) for producing a massless parton of type a in the same reaction through the 
following relation: 

^(z, Q, e) = J2 ^(*> Q> M) ® r ba (z, M , e). (2) 

b 

Eq. ((21 makes explicit the factorization of the collinear singularities present in the massless cross-section da a /dz(z, Q,e). 
All IR divergences are regularized by working in d = 4 — 2e dimensions. Throughout this paper we will work in the 
MS scheme where the collinear counterterms T a b take the form: 



2e ab 



(i) 



(3) 



Above, a s — a s (p) stands for the strong coupling constant renormalized in the MS scheme. Its relation to the 
bare couplings can be found in [3j. The renormalized coupling evolves with rif fermion flavors (including Q), /3q = 
(UC A - 4T R n f )/6 is the 0(a 2 s ) coefficient of the QCD /3-function, C A = 3, T R = 1/2 and C F = 4/3 are the QCD 
color factors and pj"" 1 ' arc the time-like splitting functions 18] . Our notations for the splitting functions follow 
Ref.[I3. 

The functions D u /q(x, fJ>/m) in Eq.([Q) are the perturbative fragmentation functions pj. They are solutions of the 
DGLAP evolution equation and satisfy an initial condition given at some scale /i = fio : 



D a/Q U, — ) 



D 1 ' 



z, ■ 



fa 



(4) 



Once the factorization scheme for the collinear singularities has been chosen, the initial condition D™ 1 is an unique, 
process independent function that has to be obtained from an additional calculation. A prescription for the evaluation 
of the initial condition can be extracted from the PFF method. The key observation is that if the initial scale fio is 
chosen to be of the order of the mass m, then the initial condition can not contain large logarithms and therefore can 
be derived from fixed order perturbative calculations. In particular, _D™ can be cast in the following form: 



V TO / ^— ' 



a s {fo) 
2tt 



TO / 



(5) 



Note that from the above equations and from the LO result d^ ~ fi a Q it also follows that the knowledge of D ml to a 
given perturbative order allows one to perform MS scheme subtraction (at the same order) in calculations where the 
collinear divergences are regulated with non-zero quark mass (see for example Ref.|l9|| for more details on that point). 

To evaluate the initial condition one needs to combine Eqs.QJ and (0) and truncate Eq.QJ to a fixed perturbative 
order by identifying the fragmentation function with its initial condition. As a result, D™ 1 can be expressed as a 
convolution of the following perturbative factors: a known universal counterterm and two cross-sections daQ and 
da a that can be separately evaluated in perturbation theory. Using results for heavy and massless quark production 
in e + e~, the functions da were first obtained within this approach 0. 
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Although this method works to any perturbative order, its application beyond next-to-leading (NLO) order is 
impractical due to the increased complexity in the evaluation of massive cross-sections at higher orders in perturbation 
theory. A better, process independent method for the calculation of the initial condition da^ was suggested in [TllEfil ] 
and further developed in Q where the functions a = Q,Q,q,q were calculated. The method is based on the 

observation that one needs to compute only a sufficient , process independent part of the cross-sections for massive 
(respectively massless) quark production. Moreover, this calculation can be formulated in a process independent way. 

(2) 

The main purpose of this paper is to evaluate the function d g using this approach and thus to complete the derivation 
of the initial condition of the PFF through order 0(a 2 ). With this result and after the three- loop time- like splitting 
functions become available, one will be able to evaluate the perturbative fragmentation function with next-to-next-to 
leading logarithmic accuracy. It follows from Eq.0J that the knowledge of the PFF with such accuracy will make it 
possible to evaluate heavy quark spectra at NNLO in any process and up to power corrections 0(m/Q) from pure 
massless calculations and resum large quasi-collinear logs In (m 2 /Q 2 ) with NNLL accuracy to all orders in the strong 
coupling constant. 

Next, we describe the process independent derivation of the gluon initiated component d*f\ As was shown in 0, 
the initial condition for the perturbative fragmentation function can be written as: 



D T (*' 3 = E T ^ <*) ® D b/Q (z, 3 , (6) 



with r a fc given in Eq.Q. 

The function D b /Q can be interpreted as a bare fragmentation function for production of a massive quark Q in the 
decay of an (off-shell) parton of type b. That function embodies all the effect of collinear radiation with q± ~ m < uq 
in that process. In general, the collinear limit is defined as the kinematic limit where the relative transverse momentum 
of two or more particles vanishes. When masses are introduced, the collinear limit is plf : 

q± = KQj_, m = nm, K — > 0. (7) 

The following Sudakov parametrization is particularly convenient for the discussion of the collinear limit in the 
process g* — > Q + X at order 0(a 2 ). Consider the one-to-three collinear splitting of a gluon. We denote the four- 
momenta of the produced massive quark and antiquark by (71,2 while q 3 denotes the momentum of the radiated gluon. 
The momentum of the decaying off-shell gluon is denoted as p : 

p = qi + q 2 + q 3 . (8) 

We parameterize the collinear direction by p and introduce another light- like vector n as the complimentary light-cone 
vector. We then write: 

n 

qi = Zip + fii- — z+qi,±- (9) 
(pn) 

The components /3i are found from the on-shell conditions q 2 = m 2 : 

l%= ' ' , i = 1,2,3, 10 

2z l 

with p 2 — 0. We are interested in the symmetric case: 

zi + z 2 + z 3 = l, qi tX + qi,i. + q3,x = 0. (11) 
Using the above equations, we express the momentum p through p and n: 

L Ti 

^ + 2*M' (12) 



where: 



02 = 2(pg 2 ) + 2(pg 3 ) - 2( g2 g 3 ) 



We often write z instead of Z\ to denote the measured energy fraction of the heavy quark in the final sate. Eas. l|9lll|) 
imply the energy conservation condition: 

z = l- ( nq2> - (" g3 ) (13) 
(pn) (pn) 
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The kinematics of the one-to-two splitting can be obtained by setting the momentum q 3 to zero in Eas. (|8ll3() . 

The explicit expression for the bare fragmentation function appearing in Eq. (jSJ follows directly from the factorization 
of phase-space and matrix elements in the collinear kinematics. Consider a hard scattering process characterized by 
some scale Q 3> m. Suppose that n + 2 partons with momenta k\,... , k n —i, 91, 92, 93 are produced. We consider 
momenta k%, . . . , fc n -i as non-exceptional, while momenta 91, 92, 93 are collinear, as described above. We denote the 
phase-space of n + 2 partons as dPS^" +2 - ) (fci, . . . , k n —i, 91, 92, 93)- In the limit when 91, q 2 and 93 become collinear and 
are resulting from the decay of an off-shell gluon, 91+92 + 93 = V + 0(q±), the (n + 2)-particle phase space factorizes: 

dPS (n+2) (fc!, . . . , 91, 92, 93) = dPS (n) (£i, . . . , fc n _i,p) ~[dq 2 ; m][dq 3 ; 0]. (14) 
The factor [dq;m q ] in Eq. l|14|l is the ci-dimensional one-particle phase space: 



[dq\m q ] = (2 ^l^ + (9 2 -^)- 



(15) 



We next discuss the factorization of matrix elements in the collinear kinematical configuration described above. 
As follows from Eq.Q, we are not interested in power-suppressed contributions to the cross-section. To identify 
the logarithmically enhanced terms we apply well known power counting arguments Q in the collinear limit for 
the scattering amplitudes: we rescale all momenta and masses according to Eq.Q, and then identify the leading 
contributions of the amplitudes in the limit k — > 0. Throughout the paper we work in the light-cone 1 gauge 
n M # = 0, n 2 — 0. This is necessary for the process-independent derivation of D im because in such gauges non- 
vanishing contributions to the fragmentation function are produced only by diagrams where collinear radiation is 
both emitted and absorbed by the same parton. In this paper we are concerned with the case of collinear decay of a 
gluon. Examples of relevant diagrams are shown on Fig.(l). 

For the derivation of the fragmentation function we need to consider only the spin-averaged case 2 . Then, in the 
collinear limit Eq.|7|), matrix element factorization holds in the following form: 

|M("+ 2 )(A; 1) ... ! fc n _ 1 ,9i,9 2 ,93)| 2 = \M^(k u . . . , k n . x ,p)\ 2 W(n,p,q 2 ,q 3 ) + O(k). (16) 

Here \M^ (...., p)\ 2 is the squared matrix element for producing n on-shell particles with non-exceptional momenta. 
All the contribution to the fragmentation function from the collinear splitting of a gluon is contained in the function 
W. The form of that function in presence of massless fermions is given in |22l| : 

The factor V^„ represents the squared, color averaged sum of diagrams contributing to the decay of an off-shell 
gluon with momentum given by Eq.JSJ). 

It is easy to generalize this result to the case when the final state quark and antiquark have non-zero mass m. 
Repeating the considerations in [22j one notices that the inclusion of non-zero masses does not introduce new tensor 
structure in the tensor decomposition of the function V^„. Then, since the limit m — > is regular, the only effect of 
the inclusion of the mass m is to alter the corresponding scalar form-factors in that tensor decomposition. Finally, 
we recall that in the collinear limit the mass m is considered to be of the order of the transverse momenta of the 
collinear particles. Therefore, we conclude that the inclusion of masses does not change the dimensionality of the 
scalar form-factors or their scaling in the collinear limit, i.e. in the presence of non-zero quark masses the function 
W is also given by Ea. (|T7|) . 

Combining Eas. (|lll3ll4ll(jfl we find that the function D 9 /q(z) appearing in Eq.© can be written as: 

D g /a(z) = \j W [dq 2 ;m][dq 3 ;0}8 (l - z - ^ - ^) . (18) 

Similar expression holds for the one-loop virtual contribution to that function. In this case one should replace in 
Eq. (|18JI the phase-space factor [dq 3 ; 0] with the corresponding phase-space for the virtual integration and omit the 
term with q 3 in the energy-conservation delta- function. In both cases the function W is evaluated from Eq. l|17|l and 
the tensor V^„ is constructed from the appropriate diagrams. Two examples of such diagrams are shown on Fig.(l). 



1 For simplicity, the gauge fixing vector n is chosen to be the same vector that appears in the Sudakov parametrization Eq.IqI. 

2 In general, in the case of decaying gluon the matrix elements do not factorize due to the presence of non-trivial spin correlations L'ill-' fl . 



Q 



Q 



FIG. 1: Examples of diagrams for the gluon decay process g — > Q + X at 0{a 2 3 ) that contribute to the function W . The dashed 
vertical line indicates the intermediate state that has to be considered. 

III. RESULTS 

The first step in the evaluation of D„/o(z) is to obtain the splitting function W in Ea. (|17|l . Our calculation 
generalizes the tree- level massless results in j^, HjjJ ■ The next step is the integration over the corresponding real 
and/or virtual momenta. To perform both the phase-space and loop integrals, we follow the methods described in 
pil I2M l2t|: we first present all non-trivial phase-space integrals as loop integrals and then use standard multiloop 
methods |22|, such as integration- by-parts and recurrence relations, to reduce all the phase-space and loop integrals 
that have to be evaluated to a set of 18 master integrals. For the reduction to master integrals we use the algorithm 
p8| implemented in |2^| . All algebraic manipulations have been performed using Maple [3(| and Form [3l| . 

The coupling renormalization requires the result for the function D 9 /q(z) at order a s including terms of order 0(e). 
The corresponding result can be found in 0. Field renormalization is performed by calculating one-loop self-energy 
insertions in the external quark legs. The use of light-cone gauge does not lead to appearance of singularities that 
are not regulated by dimensional regularization. As a last step we perform the collinear renormalization according to 
Eqs.© and ©. 

(2) 

We decompose the result for d g (z^o/m) in terms of the color factors: 

df (z, ^) = C F T R Ff ^ + C A T R Ff* T «) + TlFp + T^F^ , (19) 
where ni = nj — 1 is the number of massless flavors. The functions F g read: 

f {c f t r ) = I _ 4z ~^ Z + 1 ln(z) + z - i + (2z 2 -22+1) ln(l - z) j>L 2 + {^{Az 2 - 6z + 3)Li 2 (z) 
4: ~ ~ 21 _1 ln 2 (z) + 6(2z 2 -2z + l)ln(l - z) ln(z) - 3(2z 2 - 2z + l)ln 2 (l - z) + ^— - ln(z) 



2 w ■ \ / \ / w , v 2 

. _ , , 10z 2 -13z + 8 2tt 2 z 2 1 „,„ o „ sf-r. (Iz-l\ (2z-l 
(4z + l)ln(l-z) + ^L-16(2z 2 -2z + l) f Li 3 [—^—) + Ll 3 ( T~ f 

-2(8z 2 - Wz + 5)Li 3 (z) - 16(2z 2 - 2z + l)Li 3 (l - z) + j - 8(2z 2 - 2z + 1) ln(z) + 8(2z 2 -2z + l) ln(l - z) 

8(l-2z) 3 l (2z-\\ ( , , s , , , , , x , N 32z 3 - 24z 2 + 12z - 1 

- >Li 2 + -12(2z 2 - 2z + 1) ln(l - z) + (8z 2 - lOz + 5) ln(z) + 
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36z 2 -34z+17, 3 , , 5(2z 2 -2z + l), 3 . , 29(2z 2 - 2z + 1) , . , , 2 . 
xLi 2 (z) ln 3 (z) + ln 3 (l - z) ^ '- ln(z) ln 2 (l - z) 

ll(2z 2 - 2z+ 1) , o. .„ . 256z 3 - 372z 2 + 156z - 29 , 2 . , 32z 3 - 46z 2 + 22z - 5 , 2 , 
+ — ~ 2 " ln 2 (^) ln(l -z) ln 2 (z) ln 2 (l - z) 

64z 3 - 87z 2 + 33z + 5 , , , . /96z 3 - 196z 2 + 255z - 48 ,„ 2 „ ^7^ 2 ^ 1 , s 

+ 3 H*) ln d - z ) + 12 ( ~ + } yJ (z) 

24z 3 - 65z 2 + 79z - 26 ,„ , „ ^1tt 2 \ , s „ ■>,„ 2 „ ,,,, , , , „ s , 

(2z 2 - 22 + 1) ln(l - z) - 2tt 2 (2z 2 - 2z + 1 In z) - In 1 - z)\ 

3 6 / 11 

2tt 2 , 1 s9 / q 9 N 7T 2 , r, . . . 152z 2 - 173z + 69 , s 

+— 1 - 2z 1 - 2z) 2 + (32z 3 - 54z 2 + 6z + 19 — + (56z 2 - 60z + 30)C 3) , (20) 

3 1 1 36 12 
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9z 2 — 2z — 14 2 I I 
f (c a t r ) = I ( 4z + ^ ^ + ( 2z 2 _ 2z + 1) ln(l - z) - — + — \ L 2 + i - 8(2z 2 - 2z + l)Lia(z) 

-2(2z 2 + 2z + l)Li 2 (-z) + (6z + 1) ln 2 (z) - 4(2z 2 -2z + l) ln(z) ln(l - z) - 2(2z 2 + 2z + 1) ln(z) ln(l + z) 

. 2 o j., , 2(2z 2 + 17z + 2) , . , 10z 2 -10z+ll M , 2tt 2 z 178z 2 - 95z + 13 
+ (2z 2 - 2z + 1) ln 2 (l -z)- -i i ln(z) + ln(l - z) - — + 

-^ji + 2(2z 2 + 2z + 1) [Si, 2 (-z) - 5i, 2 (z 2 )] + 8(2z 2 -22 + 1) (li 3 (j^J^j + Li 3 (^~T 
+ (6z 2 + 18z + 7)Li 3 (z) + 2(2z 2 - lOz + l)Li 3 (l - z) + (2z 2 + 2z + l)Li 3 (-z) + \ 4(2z 2 - 2z + 1) ln(z) 



-4(2z 2 - 2z + 1) ln(l - z) - 4 ^ 2z ^ }U 2 (— — - ) + I - (6z 2 + 2z + 5) ln(z) + 2(2z 2 - lOz + 1) ln(l - z) 



3 I V z 

-4(2z 2 + 2z + 1) ln(l + z) - — — _ + _ }LUz) + \ - (2z 2 + 2z + 1) ln(z) - 2(2z 2 + 2z + 1) 

J I 

x ln(l + z) + 2z(l + z)\u 2 (-z) + 8 ^ + 2Z + 5 ln 3 (z) - 7 ( 2z " ~ 2z + ^ ln 3 (1 _ z) + r 1Qz 2 _ 18z + 5) 
16 6 

2z 2 + 27+1 

x ln(z) ln 2 (l - z) - 3(2z 2 -2z + l) ln 2 (z) ln(l - z) ln 2 (z) ln(l + z) - (2z 2 + 2z + 1) 

i ( \ i 2/i x 32z 3 -50z 2 + 25z-6 1 2 , , /64z 3 - 44z 2 + lOz - 15 4V,.... 

x ln(z) ln 2 (l + z) H In 2 (z) ln(z) ln(l - z) 

6 \ 6 3z/ 

. . . 48z 3 - 44z 2 - 16z + 23 , 2 . . /24z 3 - 624z 2 - 475z - 181 
+2z(l + z) ln(z) ln(l + z) + ln 2 (l - z) - ^ 

, n 2 , , . /\ „ 2 oN 7T 2 12z 3 + 76z 2 -28z-22\ , 
+ (2z 2 + 6z + 1)— J ln(z) + I (22z - 6z 2 - 3)— + J ln(l - z) 

2 2 

+y (2z 2 + 2z + 1) ln(l + z) + tt 2 (2z 2 - 2z + 1)| ln(z) - ln(l - z) \ - y |l - 2z\{\ - 2zf 

-<!*■ + 24^ + 12, + »)£ - 2(10^ - 3= + 7)«3) - g^iggi^Hl + » , (21) 

36 54 27z 



(2 -) 2(2z 2 -2z + l) r2 f 4(2z 2 -2z+l) n 4(4z 2 - 4z + 5) t f 
^9 = 3 ^ + <j 3 OOO + M 1 - «)) g } L 

2z 2 -2z+l 



(ln 2 (z) + ln 2 (l - z)) + 2(2z 2 - 2z + 1) ln(z) ln(l - z) 



3 

128z 5 - 320z 4 + 160z 3 + 40z 2 - 160z + 50 , s 128z 5 - 320z 4 + 160z 3 + 120z 2 - 240z + 102 , 
45 Hz) + ^ ln(l - z) 

,„ , „ . 7T 2 384z 4 - 768z 3 + 568z 2 - 184z - 280 

-(2z 2 -2z + l)— + — , (22) 



F?^ = - 2(2Z2 ~ 2Z+1) X 2 + (- 4 ^- 2g + 1 V W + ln(l - z)) - 4 ^ 2 - 42 + 5 H L - 2Z2 ~ 2Z+1 



_ , , . Nl2 2(4z 2 - 4z + 5) r , , N , , sl , , . 7r 2 88z 2 - 88z + 56 
x [ln(z) + ln(l - z)] - — [ln(z) + ln(l - z)] + (2z 2 - 2z + 1)- — , (23) 



3 27 



where L = \a{^/ra 2 ). 

Eg. i|T$|) is the main result of this paper. 
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We next comment on the properties of Ea. (|19fl . It contains polylogarithmic functions up to rank three. These 
functions are defined through: 

Li n (z) = [' Li "- l(:r) dx, Lix(z) = - ln(l - z), Si, 2 (*) = ~ [' ^ ~ x) dx. (24) 
Jo x 2 J x 

The mass-enhanced terms proportional to ln fc (/j,§/m 2 ) , k = 1, 2 are of the general form discussed in The behavior 
of Eq. ll9|) in the limits z — ► and z — > 1 is as follows: in the limit z — > 1 the gluon initiated initial condition is only 
as singular as ~ ln 3 (l — z) and therefore suppressed with an inverse power of N for large values of the Mellin variable. 
For small z, the function d g exhibits ~ CaTr/z pole and, in addition, has less singular terms ~ ln"(z) , n < 3. 

There are two more terms in d g that deserve comment. They contain the functions | ln(z) — ln(l — z)\ and 
1 — 2z\. Those terms are proportional to the color factor (Cf ~ Ca/2)Tr and originate from the virtual diagram 
representing the quark-gluon vertex correction shown on Fig.(l). To interpret these terms we first rewrite them as 
2arctanh(y / (1 — 2z) 2 ) and — 2z) 2 respectively. We see that in the point z — 1/2 these functions are not smooth. 

At the point z = 1/2 the virtual diagram in question represents a configuration where a quark-antiquark pair is 
created and the both final-state particles have the same energy. Since their relative transverse momentum is small, 
that configuration represents effectively a threshold production of a heavy pair [33|. It is well known (see [3i| for 
example) that the one-loop vertex correction for heavy pair creation close to threshold contains powers /3 n , n > — 1 
of the velocity: 



f-t-TF' (25) 



where m is the mass of the quark and Q the invariant mass of the pair. 

One can easily verify that in our collinear kinematics the analogue of (3 takes the form: 



42(1 (26) 



1 + \q 2 ± \/m 2 



and for vanishing relative transverse momentum q± between the two final-state particles the threshold condition is 
indeed z = 1/2. With an explicit calculation one can also check that terms from the virtual diagram that are of the 
form 7 A1 /3 2 ™ -1 , n > 0, indeed lead to contributions to d g 2 ^ that are proportional to the functions arctanh(y / (1 — 2z) 2 ) 
and - 2z) 2 . 



IV. CONCLUSIONS 



In this paper we compute the gluon initiated contribution to the initial condition for the perturbative fragmentation 
function of a heavy quark through order 0(a 2 ) in the MS scheme. To derive this result we follow the general approach 
developed in The result has the following properties: the terms with mass logarithms are of the general form 
predicted by the DGLAP equation. We find that the gluon initiated contribution is suppressed with an inverse power 
of N in the limit N — > oo, where N is the conjugated to z Mellin variable. The result also exhibits a non-smooth 
behavior at the mid-point z = 1/2 which we interpret as due to Coulomb interaction at threshold. Unlike the order 
0{a s ) result, the gluon initiated initial condition at order 0(a 2 ) does not vanish when m = fiQ. Our calculation 
shows the presence of non-trivial constant term. 

Combined with the results for the fermion-initiated contributions derived in [3j this result completes the calculation 
of the initial condition for the perturbative fragmentation function through order 0(a 2 ). After the three-loop time- like 
splitting functions become available, our results will permit the resummation with next-to-next-to-leading logarithmic 
(NNLL) accuracy of the large quasi-collinear logs ~ ln fc (m 2 ) that appear in heavy quark differential distributions, i.e. 
will extend the Perturbative Fragmentation Function formalism 0] to the NNLL level. 

Our result has large number of potential applications due to the process independence of the PFF method. Of 
central importance will be its application to e + e~ in order to extract with NNLL accuracy the non-perturbative 
fragmentation components. With such result at hand, one can make high precision predictions for heavy flavor energy 
and transverse momentum distributions at hadron colliders, top-decay, DIS, etc. Combined with the previously 
calculated fermion initiated components, this result can also be used to evaluate spectra of heavy quarks at NNLO 
in any process and up to power corrections (m/Q) n , from purely massless calculations. 
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